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In the approach of the effective field theory of modified gravity, we derive the equations of motion 
for linear perturbations in the presence of a barotropic perfect fluid on the flat isotropic cosmological 
background. In a simple version of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, which is the 
minimum extension of Horndeski theories, we show that a slight deviation of the tensor propagation 
speed squared c? from 1 generally leads to the large modification to the propagation speed squared 
c 2 of a scalar degree of freedom </>. This problem persists whenever the kinetic energy px of the 
field f> is much smaller than the background energy density p m , which is the case for most of dark 
energy models in the asymptotic past. Since the scaling solution characterized by the constant ratio 
px/pm is one way out for avoiding such a problem, we study the evolution of perturbations for a 
scaling dark energy model in the framework of GLPV theories in the Jordan frame. Provided the 
oscillating mode of scalar perturbations is fine-tuned so that it is initially suppressed, the anisotropic 
parameter tj = —<3?between the two gravitational potentials Hi and significantly deviates from 
1 for cf away from 1. For other general initial conditions, the deviation of cf from 1 gives rise to the 
large oscillation of Hi with the frequency related to c 2 . In both cases, the model can leave distinct 
imprints for the observations of CMB and weak lensing. 


I. INTRODUCTION 

The constantly accumulating observational evidence for the late-time acceleration of the Universe M implies that 
there may be at least one additional degree of freedom to the system of the Einstein-Hilbert action with non-relativistic 
matter and radiation. One simple example is a canonical scalar field 4> with a sufficiently flat potential V (</>)- dubbed 
quintessence Q. The A-Cold-Dark-Matter (ACDM) model corresponds to the non-propagating limit of quintessence 
(i.e., the vanishing kinetic energy) with the constant potential A. 

The scalar degree of freedom also arises in modified gravitational theories as a result of the breaking of gauge 
symmetries present in General Relativity (GR) [gj. In f(R) gravity, for example, the presence of non-linear terms in 
the 4-dimensional Ricci scalar R gives rise to the propagation of an extra gravitational scalar degree of freedom dubbed 
scalarons HI- Provided that the functional form of f(R) is well designed it is possible to realize the late-time 
cosmic acceleration, while suppressing the propagation of the fifth force in regions of the high density through the 
chameleon mechanism El- „ 

Another well-known example of single-scalar modified gravitational theories is the covariant Galileon | [l2| , in which 
the field derivatives have couplings with the Ricci scalar R and the Einstein tensor G^ v (see Ref. [l3| for the original 
Minkowski Galileon). In this case the field kinetic terms drive the cosmic acceleration q m, while recovering the 
General Relativistic behavior in local regions through the Vainshtein mechanism [l6( . 

Many dark energy models proposed in the literature (including quintessence, f(R) gravity, and covariant Galileons) 
can be accommodated in Horndeski theories 0 -most general second-order scalar-tensor theories with a single scalar 
field (j) (see also Refs. 0 )- Using the linear perturbation equations of motion on the flat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) background [l9|, the dark energy models in the framework of Horndeski theories can be 
confronted with the observations of large-scale structure, CMB, and weak lensing [2Q 1. In such general theories the 
screening mechanisms of the fifth force in local regions were also studied in Refs. [2l| . 

There is another approach to the unified description of modified gravity- the effective field theory (EFT) of cosmo¬ 
logical perturbations [HI- [HI]. By means of the Arnowitt-Deser-Misner (ADM) formalism with the 3+1 decomposition 
of space-time [37}, one can construct several geometric scalar quantities from the extrinsic curvature K^ and the 
3-dimensional intrinsic curvature IZ^, e.g., K = A'%, S = K /IV K flv , 7Z = TZ^^. The EFT of modified gravity is based 
on the expansion of a general Lagrangian L in unitary gauge that depends on these geometric scalars, the lapse N, 
and the time t. In fact, Horndeski theories can be encompassed in such a general framework with several conditions 
imposed for the elimination of spatial derivatives higher than second order [HI ■ 

Recently, Gleyzes et al. 0] proposed a generalized version of Horndeski theories by extending the Horndeski 
Lagrangian in the ADM form such that two additional constraints are not imposed. For example the Horndeski 
Lagrangian involves the term L 4 = A^K 2 — S) + BfTZ, where the functions A 4 and B 4 , which depend on (f> and 
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X = d^cf>d M </>, have a particular relation A 4 = 2X(dBi/dX) — B 4 . In GR the ADM decomposition of the Einstein- 
Hilbert term (M 2 l /2)R, where M v \ is the reduced Planck mass, leads to the Lagrangian L 4 with B 4 = —A 4 = M 2 l /2. 
The theories with B 4 ^ — A 4 belong to a class of GLPV theories. 

On the flat FLRW background, the Hamiltonian analysis based on linear cosmological perturbations shows that 
GLPV theories have only one scalar propagating degree of freedom [35|, • One distinguished feature of GLPV 

theories is that the scalar and matter sound speeds are coupled to each other (.'i ll. [35|. For example, in the covariantized 
version of the Minkowski Galileon where partial derivatives in the Lagrangian are replaced by covariant derivatives, 
the scalar propagation speed squared Cg becomes negative in the matter-dominated epoch due to a non-trivial kinetic- 
type coupling [411 ]. Unlike the covariant Galileon with positive c 2 during the matter era, the covariantized Galileon 
mentioned above (a class of GLPV theories) is practically excluded as a viable dark energy scenario. 

In this paper, we develop the analysis of cosmological perturbations further to confront dark energy models in 
GLPV theories with observations. First, the linear perturbation equations of motion are derived in the presence of 
a barotropic perfect fluid for a general Lagrangian encompassing GLPV theories. We provide a convenient analytic 
formula for c 2 and show that even a slight deviation from Horndeski theories generally gives rise to a non-negligible 
modification to the scalar sound speed. 

We also apply our general formalism to a simple dark energy model with a canonical scalar field (f> in which the 
function B 4 differs from — A 4 = M 2 x /2. In this case the tensor propagation speed squared c 2 = —B 4 /A 4 is different 
from 1. We show that this deviation leads to a significant modification to Cg whenever the field kinetic energy px is 
suppressed relative to the background energy density p m . The scaling solution characterized by the constant px/Pm 
is a possible way out to avoid having large values of in the early cosmological epoch. 

For the scaling dark energy model described by the potential V(<f>) = V\e~ Xl ^^ Mpl + I/ 2 e -A2< ^ Mpl (Ai > 10 and 
A 2 1), we study the evolution of cosmological perturbations and resulting observational consequences. For the 
initial conditions where the contribution of the oscillating mode Vm ' 1 to the velocity potential V m is suppressed, the 
evolution of perturbations is analytically known during the scaling matter era. In particular, the anisotropic parameter 
77 = — <f > /'F between the two gravitational potentials and $ exhibits a large deviation from 1 for c 2 away from 1. 
If the oscillating mode Vm 1 gives a non-negligible contribution to V m initially, the rapid oscillations with frequencies 
related to c 2 arise for the perturbations like ’F and V m . Thus the model in the framework of GLPV theories can be 
clearly distinguished from that in Horndeski theories. 

This paper is organized as follows. In Sec.|TT]the extension of Horndeski theories to GLPV theories is briefly reviewed. 
In Sec. Illll the perturbation equations of motion are derived in the presence of a barotropic perfect fluid according to 
the EFT approach encompassing GLPV theories. In Sec. HVI we obtain convenient formulae for the matter and scalar 
propagation speeds in the small-scale limit. In Sec. |V] we present observables associated with the measurements of 
large-scale structure, CMB, weak lensing, and discuss the quasi-static approximation on sub-horizon scales. In Sec.lVD 
we propose a simple dark energy model in the framework of GLPV theories and study its observational signatures by 
carefully paying attention to the oscillating mode induced by Cg. Sec. IVIIl is devoted to conclusions. 


II. THE THEORIES BEYOND HORNDESKI 

The EFT of cosmological perturbations is a powerful framework to deal with low-energy degrees of freedom in a sys¬ 
tematic and unified way for a wide variety of modified gravity theories. It is based upon the 3+1 ADM decomposition 
of space-time described by the line element [37| 

ds 2 = g llu dx tJ ‘dx 1 ' = —N 2 dt 2 + hij{dx l + N l dt){dx - 7 + N^dt) , (2.1) 

where N is the lapse function, N l is the shift vector, and hij is the three-dimensional spatial metric. The extrinsic 
curvature is defined by K = /i A n„ ; A, where = (— N, 0,0, 0) is a unit vector orthogonal to the constant t hyper¬ 
surfaces Ej and a semicolon represents a covariant derivative. We also introduce the three-dimensional Ricci tensor 
TZfiu = ^Rp V on £ t . Then we can construct a number of geometric scalar quantities: 

I< = a % , s = K^K^, n = n%, z = n, w n>‘ v , u = u^k^. ( 2 . 2 ) 

Horndeski theories [f3] are the most general scalar-tensor theories with second-order equations of motion in generic 
space-time. The action of Horndeski theories is given by S = f d 4 x^/—g L with the Lagrangian El 

L = G 2 {ct>,X) + G 3 ^,X)n<t) + G4(.^X)R-2G4,x(.&X) [(O0) 2 - 

+G^X)G I1 ^ + iG 5 ,x(</>,*)[(D^) 3 -3(0^)^^ + 2$.^"ft".„] , 


(2.3) 
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where □</> = and the four functions G, (* = 2,3,4, 5) depend on a scalar field <f> and its kinetic energy 

X = g^d^d^cj), with G it x = dGi/dX. Choosing the unitary gauge (f> = </>(t) on the flat FLRW background described 
by the line element ds 2 = — dt 2 + a 2 (t)Sijdx l dx ^, the Lagrangian (12.311 can be expressed in terms of the geometric 
scalars introduced above, as [29[ 

L = A 2 {N 7 1) + A 3 (N, t)K + A 4 (N, t)(K 2 - S) + B 4 {N, t)K + A 5 {N , t)K 3 + B 5 (N, t) (U - K1Z/2) , (2.4) 

where K 3 = K 3 — 3KK flv K 111 ' + 2K^ v K fJ,x K v Up to quadratic order in the perturbations we have K 3 = 3H(2H 2 — 
2KH + K 2 — S), where H = d/a is the Hubble parameter (a dot represents a derivative with respect to t). 
Horndeski theories satisfy the following two conditions [29 ] 

A 4 = 2 XB 4 , x ~B 4 , A 5 = -±XB 5 , x . (2.5) 

More concretely, the coefficients Gi in Eq. (12.31) and A tl B, in Eq. (12.41) are related with each other, as 

A 2 = G 2 - XF 3 ^ , A 3 = 2(—X) 3 / 2 F 3i x ~ 2 yf^XG^ , 

A 4 = —G 4 + 2XG 4t x + XGs^/ 2 , B 4 = G 4 + X(G§ t( j, — F§^/ 2 , 

A 5 = -(-Xf/ 2 G 5 , x /3, B 5 = -V=XF S , (2.6) 

where F 3 and F 5 are auxiliary functions obeying the relations G 3 = F 3 + 2XF 3 x and G 3 .x = F$/(2X) + F 3 _x- Since 
X = —cf> 2 (t)/N 2 in unitary gauge, the functional dependence of <f> and X can translate to that of t and N. 

It is possible to go beyond the Horndeski domain without imposing the two conditions (I2.5[) [35) . This generally 
gives rise to derivatives higher than second order, but it does not necessarily mean that an extra propagating degree of 
freedom is present. In fact, the Hamiltonian analysis on the flat FLRW background shows that the theories described 
by the Lagrangian m , dubbed GLPV theories, do not possess an extra scalar mode of the propagation (35i. l38l - [40| . 


III. PERTURBATION EQUATIONS OF MOTION 

The Lagrangian (12.41) depends on N, K, S, 1Z,U,t, but not on Z. The dependence on Z ap pears in the theories 
with spatial derivatives higher than second order [33j, [13, e.g., in Horava-Lifshitz gravity (42j| . In the following we 
shall focus on the theories described by the action 

S = J d A Xy/^gL(N,K,S,n,U;t) + J d^xyf^g , (3.1) 

where L m is the Lagrangian of the matter field 'F m . We consider a metric frame in which the scalar field (j) is 
not directly coupled to matter (dubbed the Jordan frame). For the matter component we consider a scalar field \ 
characterized by 

L m = P(Y), Y = g^d f ,. X d,X, (3.2) 

whose description is the same as that of the barotropic perfect fluid (43I - I45) . The perfect fluids of radiation and non- 
relativistic matter can be modeled by P{Y) = ciF 2 and P(Y) = c 2 (Y — Y 0 ) 2 with |(Y — lo)/Yo| ^ respectively, 
where ci,C 2 ,T 0 are constants I4lll46). 

The linearly perturbed line element on the flat FLRW background with four metric perturbations A , ip, </, E and 


•J ' ' t J J 

tensor perturbations 7 ^ is given by [43] 

ds 2 = — (1 + 2 A)dt 2 + 2 d$dtdx' 1 + a 2 (t) [(1 + 2 ()6ij + 28$jE + 7 ^] dx l dx 2 . (3.3) 

Then the shift vector Ay is related to the perturbation t/>, as 

Ni = dii >. (3.4) 

In the following we choose the gauge conditions 

5(j) = 0, E = 0, (3.5) 

under which temporal and spatial components of the gauge-transformation vector are fixed. The background values 
(denoted by an overbar) of the ADM geometric quantities are 

= Hh^ , K = 3H , S = 3H 2 , n, w = 0, H=U = 0. (3.6) 
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Around the FLRW background we consider the scalar perturbations SN = N— 1, SK^ = K^—Hh^, SS = 2H5K+ 
5K^^8K V and 1Z = 8 {R + 5 2 R, where 5iTZ and S 2 1Z are the first-order and second-order perturbations respectively. 
The scalar U. which is a perturbed quantity itself, obeys the relation f d 4 XyJ—ga(t)U = f d 4 Xy/—g[a(t)lZK/2 + 
a(t)lZ/(2N)\ up to a boundary term, where a(t) is an arbitrary function with respect to t. 

Decomposing the scalar field \ as X = X.(t) + $x(t> x ) and omitting the overbar in the following discussion, the 
kinetic term Y, expanded up to second order, can be expressed in the form Y = — \ 2 + <5i Y + S 2 Y , where 

SiY = 2 x 2 SN-2x6 X , (3.7) 

feV = -5 X 2 -3ic 2 SN 2 +4x6 X 6N + % 6 ^d^d J 6 x+ -^(dS X ) 2 , (3.8) 

a z a z 

and (cMx ) 2 = S ' 13 d,SxdjSx- The energy-momentum tensor of the field x is given by = Pg^ u — 2 
Defining the linear perturbations of energy density, momentum, and pressure, respectively, as ST® = — Sp , STf = diSq , 
and 5Tj = 8 PS*, it follows that 

Sp = {P, Y + 2 YP :Y y) 6 {Y , 6 q = 2P iY x6x, SP = P, y SiY . (3.9) 

These quantities appear in the perturbation equations of motion presented later. 


A. Background equations of motion 


Ex pan ding the action m up to linear order in scalar perturbations, we obtain the first-order action = f d 4 x C\ 
with |29l . l34j 


£, = a 3 (L + P N - 3 HP - p) SN + 3 (l - T - 3 HP + p) a 2 Sa - 2 a 3 P x xSx + a^SS^R , (3.10) 

where 

T = L, k + 2HL,s, (3.11) 

£ = Ljz + —Ljj + ^HL u , (3-12) 

p = 2 YPy-P. (3.13) 

The last term of Eq. (13.101) is a total derivative irrelevant to the background dynamics. Varying Eq. (13.101) with 
respect to SN, 8a, and Sx, respectively, we obtain the background equations of motion 

L + L jN ~ 3HT= p, (3.14) 

L- T -3HP = -P , (3.15) 

^(« 3 Pvx)= 0. (3.16) 

These correspond to the Hamiltonian constraint, the momentum constraint, and the equation of motion for Xi re_ 
spectively. Equation (13.161) is equivalent to the continuity equation p + 3 H(p + P) = 0 by using the definition (13.131) 
of the field energy density. 


B. Perturbation equations of motion 

We expand the action m up to second order in scalar perturbations to derive the linear perturbation equations 
of motion. In doing so, we use the following properties 

SK) = (C - HSN) S] - ^5 lk (< 9 k Nj + djN k ) , (3.17) 

SKij = -(Sij&C + didj C) , (3.18) 

where d 2 C, = 8 l 3 didjC,. Expansion of the action (13.11) gives rise to the terms in the forms (d 2 ip) 2 /a 4 , (d 2 , ip)(d 2 Q/a 4 , 
and ( d 2 C,) 2 /a 4 , which generate the spatial derivatives higher than second order. The absence of these higher-order 
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terms requires that 


L,kk + ^HL^sk + 4H 2 L,ss + ZL'S = 0 , 

L.kiz + 2 HL t sn + - L'U + HL^ku + 2 H 2 L t su = 0 , 


'KTl 


+ 2HL 


,KU 


H L^uu = 0 . 


(3.19) 

(3.20) 

(3.21) 


The GLPV Lagrangian (12.41) obeys all these conditions, so we assume the conditions (I3.19I) - (I3.21I) in the following 
discussion. 

The second-order Lagrangian density of C m = yf^gP{Y) = Ny / hP(Y) reads 

4? = —pSNSVh - 2 P )YX S X SVh + a 3 [p )Y S 2 Y + i P xy 5iY 2 + P^SNSiY^J , (3.22) 

where SVh = 3 a 3 C,. The Lagrangian density C = y/—gL contains the second-order contribution ( L + L >x — 
3HP)SN5\/h. On using the background Eq. (13.141) . this cancels the first term of Eq. (|3.22l) . Then, expansion of 
the total action (13.11) leads to the second-order action = f d 4 x£ 4 with 


4? = a 3 


L, n + -L nn - 3 H{W - 2 L S H) \ SN 2 + W 3(-f - 4(X> + £)-± \ SN + 4 L s 


d 2 i) 


d 2 ( 


id 2 ip 


M) 


d 2 il> 


where 


-6 L,sC 2 + - ( P V + 2 YP xy ){ X SN - 5 X Y - QP X X^xC~ 2 P,yxSx~t + P T 


W = L t KN + 2HL'SN + 4:L,sHi 
P = L,ntz — t;L,u + HL^nu . 


(dSxf 


, (3.23) 

(3.24) 

(3.25) 


Varying Eq. (13.231) with respect to SN, d 2 i/j, (, and $X respectively, and employing Eqs. (13.91) and (13.16|l . we obtain 
the following linear perturbation equations of motion 


(2L jv + L, nn - 6 HW + 12L X H 2 ) SN + - ^ 4 ) W - 4 (V + £)^-=8p, 

WSN -4L, s ( = -6q, 

44 (« 3 ^) + 4 ( P + e )^T~ + 44 2 C + 6 P,YX 2 SN = 3 SP , 

CL CLl CL CL 

44 P( p v + 2 YP, yy )(6x - X<W)1 + 3P >y xC ~ xP,y^~ - Py^ = 0, 

a 6 at L J or or 


where 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 


We note that Eqs. (13.311) and (13.321) also follow from the continuity equations ST^i-^ = 0 and ST^ o ;/i = 0, respectively. 
Substituting Eq. (13.301) into Eq. (13.281) and using Eq. (13.311) . it follows that 


d 2 ib 

y = 4L 5 — j- — 3Sq. 

a z 

On using Eq. (13.161) . it is easy to show that the momentum perturbation 5q = 2 P, yX 8 x obeys 

5q + 3H5q = -(p+ P)SN - SP. 

Similarly, Eq. (13.291) can be expressed in the following form 

Sp + 3 H (Sp + SP) = -{p + P)( 3C - ^ d ~ 6<1 


(i,s + HL s )*I> + L s ij> + (V + £)5N + £( = 0, 


(3.33) 





































6 


where the integration constant is set to 0. The dynamics of scalar perturbations is known by solving Eqs. (13.261) . 
(13.271) . (13.331) together with Eqs. (13.311) and (13.321) . 

For the tensor perturbation 7 ^ the second-order action reads = J d 4 xC^\ where [29L [ 32 L EH 


4 2> =« 




Here, the propagation speed squared is 


(3.34) 


c t = 


(3.35) 


Then the equation of motion for gravitational waves is given by 


7ii + I 3H+ L s ) Hj c t a2 


2W 1x3 = 0 . 


(3.36) 


Provided that Ly > 0 and c 2 > 0, the ghost and Laplacian instabilities are absent for the tensor mode. 


IV. PROPAGATION SPEEDS OF SCALAR PERTURBATIONS 

We derive the propagation speeds of the gravitational scalar and the matter field in the small-scale limit. In doing 
so, we first express 6 N and d 2 t/j/a 2 in terms of C, and their derivatives by using Eqs. (13.261) and (13.271) . Substituting 
these relations into Eq. (13.231) , the Lagrangian density can be expressed in the form [3J, |4l| 


(x* 


C 2 = a 6 (X f KX - djX t Gd J X - X t BX - X t MXJ , 


(4.1) 


where AT, G, B , M are 2x2 matrices and X 4 = (Q 6x/M p \). The components of the two matrices K and G are 
given, respectively, by 


K n =Q S + 


16 L' 


5 ,-,2 


M p 2 ,W 2 


X z K 22 , I<22 = {2± z P,yy - P,y)M z , 


K - W ALs * K 

A 12 — A 21 — — , r . ,J ^22 ■ 


M pl W 


G n = 2(M + HM - £), 


where 


G22 = -PvMp, : 


8 L 2 s, 


Gi 2 “ G 21 - _ l^ G22 ’ 


Qs = 6L s + 777t(2-^W + L,nn — 6 HW + \2H z Ls ), 
M = 


W 2 
4 L s (T> + £) 


W 


(4.2) 

(4.3) 

(4.4) 


The scalar ghosts are absent as long as the determinants of principal sub-matrices of K are positive, which translates 
to the two conditions Kn > 0 and Q S K 22 > 0. These conditions are satisfied for Q s > 0 and K 22 > 0. 

The dispersion relation following from Eq. dUD in the limit of the large wave number k is given by 
det ( uj 2 K — k 2 G/a 2 ) = 0. The scalar propagation speed c s , which is related to the frequency w as w 2 = c 2 k 2 /a 2 
obeys 


(cg-Kii — Gn) (c^ K22 — G22) — (c 2 K\2 — G12) = 0 . 


(4.5) 


In Horndeski theories there is the specific relation T> + £ = L,s and hence G 12 //L 12 = 022 ^ 22 - In this case the two 
solutions to c 2 are given by 


c 2 = 


Cur — 


G 


22 


P 


,Y 


SP 


K 22 Py - 2 x 2 P,yy S P 


(4.6) 


Qs 


G 11 — {Kii — Q s ) 


G 2 

K- 


22 


2 ( ■ 8 L 2 s x 2 Py' 


(4.7) 
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In GLPV theories the relation T> + £ = L s no longer holds, so we define the parameter 


qh = 


V + £ 


~ 1 = ct 1 


D 

L,s 


(4.8) 


which characterizes the deviation from Horndeski theories. We also introduce the following quantity 

;(P + P) 


Kn 

Qs 


Ph = 2 c m — -1 Uh = 


16 L 2 s y 


W 2 Q S 


-a h 


(4.9) 


If T> = 0, then the parameter an is simply related to the deviation of the tensor propagation speed squared from 1, 
as an = cj? — 1. This is the case for the theories with B 4 = B^(t) and constant B$. In Sec. IVII we shall discuss the 
dynamics of cosmological perturbations for a simple model satisfying the condition T> = 0. 

Eliminating the terms G 22 , G 11 , G 12 , and Kn in Eq. (|4.5P with the help of Eqs. (|4.6[) , (|4.7|) , (|4.9|) and the relation 
G 12 /K 12 = (1 + o;h)G 22 /^ 22 , the two solutions to Eq. (|4.5|) can be expressed as 


c 


2 

s 


c 2 m + 4 - Ph + \J {c 2 m -C 2 J + /? H ) 2 + 2 c^«h/ 3 h 
c m T C H Ph \/(i - c\ + /?h ) 2 + 2 c^a H /?H 


(4.10) 

(4.11) 


For non-relativistic matter with = 0, Eqs. (14.101) and (14.111) reduce to = 0 and c 2 = — (3 ^ respectively. 

For the general perfect fluid with ^ 0, we consider the case in which the deviation from Horndeski theories is 
small, i.e., |an| 1. Then the propagation speeds (14.101) and (14.111) are approximately given, respectively, by 


— C™ 


2 (c| - (?m - Ph) 


UhPh ) 


— c h ~ Ph + 


2(ch -c 2 m~ Ph) 


anPn ■ 


(4.12) 

(4.13) 


The condition, jap I "C 1, does not necessarily mean that |/3h| is also much smaller than 1. In the covariantized 
Galileon model [35} where the partial derivatives of the original Minkowski Galileon [l3| are replaced by the covariant 
derivatives, we have /3 h = 3/10 and c|, = 11/40 during the matter era for late-time tracking solutions (in which 
regime |an| is much smaller than 1) |4l[. The reason why |/3h| is not as small as |an| comes from the fact that the 
variable Q S /M 2 X in Eq. (14.91) is much smaller than 1 in the early cosmological epoch. Since c 2 = —1/40 in this case, 
the covariantized Galileon is plagued by the Laplacian instability on small scales. On the other hand, for the covariant 
Galileon [l8[, we have c 2 = = 1/40 during the matter era [TB}- 

Generally, the sound speed squared c 2 is subject to the modification arising from the deviation from Horndeski 
theories, such that c 2 ~ — /3h- Meanwhile, provided that |q;h/5h| "C 1, the correction to the matter sound speed 

squared c^, i.e., the second term on the r.h.s. of Eq. (14.121) . is suppressed to be small. This shows that the effect 
beyond Horndeski theories arises for the scalar sound speed c s rather than the matter sound speed c m . In Sec. IVII we 
shall apply the results in this section to a concrete model that belongs to a class of GLPV theories. 


V. CONFRONTATIONS WITH OBSERVATIONS 

In this section we discuss several physical quantities associated with the measurements of large-scale structures, 
CMB, and weak lensing in order to confront GLPV theories with observations. We then proceed to the discussion of 
the quasi-static approximation for the perturbations deep inside the sound horizon. 


A. Observables 

We first introduce the gauge-invariant combinations of the matter density contrast S m and the velocity perturbation 
Vm ; ns 


S m = 5 — 3 Hv 


Vm =V+(l+w)^- 


(5.1) 






















where 


(5.2) 


<5 = 


5p_ 

P ’ 




Since we choose the unitary gauge (5<j> = 0), the perturbation v m is equivalent to v itself. In Fourier space we can 
rewrite Eqs. (13.311) and (13.321) . respectively, as 


v m + 3 H ( c 2 m -w)v m = -(1 + w)6N - c 2 m 5 m , 

S m + 3 (Hv m y + 3 H (c 2 , - w) (6 m + 3 Hv m ) = -(1 + w) ^3C + 



(5.3) 

(5.4) 


where c m is defined by Eq. 6H). 

Since we are interested in the growth of structures after the onset of the matter era, we shall focus on the case 
of non-relativistic matter characterized by w = 0 and = 0. Taking the time derivative of Eq. (15.41) and using 
Eq. (15.31) . we obtain 

k 2 

5 m + 2H5 m + -f* = -3B-6HB, (5.5) 

or 

where B = £ + Hv m , and \F is the gauge-invariant gravitational potential defined by H3] 


\F = SN + ip . (5.6) 

If c m is not exactly 0 and the term — c^<5 m on the r.h.s. of Eq. (15.31) is not neglected, this gives rise to the term 
Cm{k 2 /a 2 )S m on the l.h.s. of Eq. (15.51) . This works as a pressure that prevents the gravitational growth induced 
by the source term (fc 2 /a 2 )\F. The matter propagation speed squared in GLPV theories is not equivalent to c^, 
but, in the limit —> 0, they are identical to each other. On the other hand, the scalar sound speed squared c 2 is 
generally subject to a non-negligible change even by the slight deviation from Horndeski theories. 

In order to know the evolution of the matter density contrast 5 m , we need to relate the gravitational potential 'F 
in (1531) with 5 m - Usually, this relation is expressed in the following form 

—~o ^ = —47T GeftpSm , (5.7) 

a z 

where G e g is the effective gravitational coupling. In GR, G e g is equivalent to the Newton gravitational constant 
G. In modified gravitational theories, G e s generally differs from G. In Horndeski theories, for example, the quasi¬ 
static approximation provides the analytic expression of G e g for the perturbations deep inside the sound horizon El- 
Provided that the terms on the r.h.s. of Eq. (15.51) are negligible compared to those on the l.h.s., the evolution of 5 m 
is known by integrating Eq. (15.51) . The growth rate of 5 m is related to the peculiar velocity of galaxies [49|. For the 
observations of redshift-space distortions, the quantity fas is usually introduced [50j . where 


f = 


H5n 


(5.8) 


and as is the amplitude of over-density at the comoving 8h _1 Mpc scale (h is the normalized Hubble constant 
Hq = 100ft. km sec -1 Mpc -1 ). 

In order to confront modified gravity models with the observations of CMB and weak lensing, we also introduce 
the following gauge-invariant gravitational potential 


$ = C + Hip. 

The effective gravitational potential associated with the deviation of light rays is given by 

4>eff = i OF - $) . 


(5.9) 


(5.10) 


Introducing the anisotropic parameter 


$ 


V = ~ 


<F ’ 


(5.11) 


Eq. (15.101) can be expressed as <F e ff = (1 + rj)'i’/2. In GR we have r) = 1 and hence <F c ff = T. We caution that the 
definition (15.111) is valid only for >P ^ 0. If the gravitational potential ^F crosses 0 with oscillations, we should compute 
<F e g from Eq. (15.101) rather than using rj. As we will see in Sec. IVI1 the crossing of = 0 can actually occur in GLPV 
theories if the oscillating mode initially dominates the perturbation v m . 
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B. The quasi-static approximation on sub-horizon scales 


For the observations of large-scale structure and weak lensing, we are primarily interested in the evolution of 
perturbations for the inodes deep inside the sound horizon ( c s k aH). In the presence of a propagating scalar 
degree of freedom, there is an oscillating mode of the field perturbation in addition to the mode induced by matter 
perturbations. Provided that the oscillating mode of perturbations is suppressed relative to the matter-induced mode, 
the time derivatives of metric perturbations (like £ and ip) can be neglected relative to the terms involving their spatial 
derivatives |52i] . In Horndeski theories, this quasi-static approximation was first employed in Ref. 0 to derive the 
analytic expression of G e g , $ e ff, and p. 

In GLPV theories, let us discuss what kind of difference from Horndeski theories arises. First of all, Eq. (13.331) can 
be written in the form 


(1 + a H )tf + cfc + 



Hip 


an ip. 


(5.12) 


Since the time derivative an ip does not vanish, we need to deal with Eq. (15.121) as the differential equation rather 
than the constraint equation. 

Under the quasi-static approximation on sub-horizon scales, the dominant contributions to Eq. (13.261) can be 
regarded as those involving the Laplacian terms d 2 ip/a 2 , d 2 (/a 2 and 5p. Then, in Fourier space, Eq. (13.261) reduces 
to 

y ^~2 ^P + 4T,s(l + an^C — P& ■ (5.13) 

ar 

Taking the time derivative of Eq. (13.261) and eliminating the term 5p on account of Eq. (13.32p . the quasi-static 
approximation for sub-horizon perturbations leads to 


VW + (W + HW + p)ip + 4 (1 + an)(l s + HLs) + L s an C --4a H i,sC- 


(5.14) 


In deriving Eq. (15.141) we have implicitly assumed that the mass of the scalar degree of freedom <p is at most of 
the order of the Hubble parameter H. In some of the modified gravity models in which the chameleon mechanism 
is at work 0, rricf, can be much larger than H as we go back to the past. In the regime H , however, the 

scalar degree of freedom is nearly frozen, so that the evolution of perturbations is similar to that in GR 10 ]. The 
modification of gravity manifests itself in the late cosmological epoch associated with the cosmic acceleration, in which 
regime is at most of the order of H. 

In Horndeski theories we have an = 0, so the terms on the r.h.s. of Eqs. (15.121) and (15.141) identically vanish. 
In this case, we can express ( and ip in terms of T by using Eqs. (|5.121) and (15.141) . Substituting these relations 
into Eq. (15.131) . we obtain the modified Poisson equation (15.71) with the effective gravitational coupling G e g. The 
effective gravitational potential 'Peg and the anisotropic parameter p are known accordingly. This procedure is given 
in Appendix A. 

In GLPV theories the time derivatives ip and ( are present, so we cannot derive the closed-form expression of G e ff 
and 4> e ff- The existence of the time-derivative terms in Eqs. (15.121) and (15.141) implies that the oscillating mode may 
play a non-trivial role for the evolution of perturbations. In Sec. IVII we shall discuss the condition under which the 
oscillating mode is suppressed relative to the matter-induced mode for a model in the framework of GLPV theories. 


VI. A CONCRETE MODEL 

In this section we consider a concrete model in which Horndeski theories are minimally extended to GLPV theories. 
The model is described by the action m with the Lagrangian 

L = A 2 + A 4 (K 2 ~ S) + B 4 ll, (6.1) 


where 

A 2 = -Ix-V{<P), A4 = —, B i = ^M^F(<p), (6.2) 

with V(cp) and F(cp) being functions of </>. The model with F(cp) = 1 correspond to GR in the presence of a minimally 
coupled scalar field <p with a potential V((p). If the function F[cp) differs from 1, then the first condition of Eq. (12.51) 
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is not satisfied, so this is already beyond the realm of Horndeski theories. For the matter Lagrangian L m we consider 
a single perfect fluid described by the constant equation of state w = P/p. 

From Eqs. lj3.14M3.15ll we obtain the background equations of motion 

3 M% i H 2 = ±j> 2 + V(<l>)+p, (6.3) 

- 2 M^H = <j ) 2 + p + P. (6.4) 

Taking the time derivative of Eq. (16.311 and using Eqs. (13.1611 and (16.411 . the scalar field </> obeys 


4>+3Hj>+Vt = 0. (6.5) 

Equations (I6.3M6.5I1 are equivalent to those in GR. This means that, at the background level, we cannot distinguish 
between the theories with same A 4 but with different B 4 flH |. 

For the later convenience we introduce the following variables [5(1 


= ^ „ = Vv 0 = P 

2:1 “ v^tfMpi ’ X 2 “ \f3HM p i ’ m ~ 3MpjiJ 2 ’ 

with which the Friedmann equation (16.31) can be written as f l m = 1 — xf — x\. 
and X 3 = 4>/M p \ obey the equations of motion 


A = — 


M p iE 0 
V 


( 6 . 6 ) 


The dimensionless quantities aq, X 2 , 


/g ^ 

x\ = -3aq + — \x\ + -aq [(1 - w)x\ + (1 + wj)(1 - xl)] , (6.7) 

x' 2 = ~^y~XxiX 2 + ^x 2 [(1 - w)xl + (l + w)(l -xl)] , (6.8) 

x '3 = \[§x \, (6-9) 

where a prime represents a derivative with respect to M = In a. Note that the variable A is known as a function of X3. 


A. Propagation speeds 


The background degeneracy between the theories with different values of B 4 is broken at the level of cosmological 
perturbations. From Eq. (13.3511 the tensor propagation speed squared is given by 




( 6 . 10 ) 


The deviation from Horndeski theories can be quantified by the difference of c 2 from 1. In fact, the deviation parameter 
qh in Eq. (HU) reads 


cxh = c t - 1. 

The quantities and /3 h defined in Eqs. (14.71) and (14.91) reduce, respectively, to 


c h — 1 + 


2 M^H 2 


2 \ ^ , 2CtCt 




, ( 2 M 2 77 n 

/?H = 2(1 - c t 2 ) 1 + — 


( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


where we used the background Eq. (10.41) . In the regime |an| 1, the sound speed squared (|4.13|) can be estimated 
as cl~ c| - /3 h , i.e., 


c 


2 

s 


l-2(l-c t 2 ) 


2 M^H 2 
(j > 2 


(1-G 2 ) 


H 
H 2 


2ctCt 

H 


(6.14) 
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The density parameters of the field kinetic energy and the potential energy are given, respectively, by fix = x\ and 
fly = In terms of these parameters, Eq. (|6~H) can be expressed as 


H 

7/2 


—3f2x 


-(1 + w)fl m . 


(6.15) 


Then, the sound speed squared (16.141) reduces to 


c 


2 

S 




(1 + w) 


2 fix 



2 c t c t 

+ 3 Hfl x ' 


(6.16) 


If the field 4> is responsible for dark energy, its kinetic energy is usually suppressed relative to the background 
energy density in the early cosmological epoch, i.e., flx/fl m "C 1 (unless we consider early dark energy models). For 
Cj > 1 the sound speed squared (16.161) can be negative, which leads to the Laplacian instability on small scales. Let us 
consider the case in which Ct = 0, i.e., F(^>) = constant. In order to realize the condition c)? > 0 during the radiation 
and early matter eras, we require that 


1<^- (6-17) 

i Lfji 

As we go back to the past the ratio flx/flm gets smaller, so c 1 needs to be very close to 1 in the early cosmological 
era. 

When cj? < 1 the Laplacian instability is absent, but c s becomes highly super-luminal. This means that the quantity 
c s k/a, which appears in the perturbation equations of motion, becomes much larger than k/a , so that the perturbation 
theory can break down for c s k/a above some cut-off scale M. 

We stress that the problem of large values of | cj | persists whenever there is a scalar field 4> whose kinetic energy is 
suppressed relative to the background energy density. The slowly varying scalar potentials V(</)) responsible for the 
late-time cosmic acceleration are generally plagued by this problem. 

One way out is to consider a primordial scaling field characterized by the constant ratio fl \/ fl m [ 5314551 ] . The scaling 
solution can be realized by the exponential potential V((j>) = y 1 e _Al,i / Mpl for the constant Ai satisfying \\ > 3(1 + w) 
[53j . For the success of the big bang nucleosynthesis, the slope Ai is constrained to be Ai > 9.4 [5g]. In this case the 
late-time cosmic acceleration is not realized, so the form of the potential needs to be modified to exit from the scaling 
regime. One of such models is given by [57] 


V(<l>) = y ie - Al ^ Mpl + y 2e - A ^/ M p‘ ) (6.18) 

where Y, V 2 , A|, A 2 are constants. Provided that A| < 2, the solutions finally approach the acceleration attractor 
characterized by the field equation of state iude = — 1 + A|/3. 


B. Perturbation equations and the oscillating mode 

Since our primary interest is the growth of structures after the onset of the matter-dominated epoch, we consider a 
non-relativistic perfect fluid characterized by w = 0 and = 0 for the matter Lagrangian L m = P(Y). For example, 
the Lagrangian P(Y) = C 2 (Y — To) 2 with |(Y — Yo)/Yo| -C 1 can describe a non-relativistic perfect fluid jUH^]. In 
the following we also focus on the case in which c £ is constant. 

First, we introduce the following dimensionless quantities 


V m = Hv m , y = Hijj. 


(6.19) 


For the model m , the perturbation equations of motion in Fourier space following from Eqs. (13.261) . (13.271) . and 














(13.31l) - (13.33l) are given by 


12 


SN + — fl m V m j 


= - 


(S 


- 1 ) X - c t 2 diV - c t 2 C , 


7 
x' 

5' = -XfI m <5 + 3(O x -l) ( W+^-^) (c^+KJ, 

c 

n x SN 


— 5N + ~jjVn J 


io A _I 
2 ^ ^mOm ^ 


fc 

7777 


( c t C + x) j 


( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 

(6.23) 

(6.24) 


where S m = 5 — 2>V m is the gauge-invariant matter perturbation. 

The two gravitational potentials m and $ are given, respectively, by 4/ = SN — (H'/H)x + x' and = C + X- From 
Eq. (16.211) we obtain the following relation 


\E' + $ = (l-c t 2 )(<W + C). 

Taking the A f derivative of Eqs. (16. 2211 - (16.231) and using other equations of motion, it follows that 


ir_ 


r) ^•m^m 


k 

Tff 


TTf 

(c t 2 - 1)SN - 12 Q X 6N + 18 n x ( 1 + — + 


n 




H 2Q X J 


V m 


(6.25) 


(6.26) 


Since = — <f> for c 2 = 1, the anisotropic parameter 77 defined by Eq. (15.111) is equivalent to 1. In this case the first 
term on the r.h.s. of Eq. (16.261) also vanishes. Provided f l x -C 1, the remaining terms on the r.h.s. of Eq. (16.261) are 
suppressed relative to those on the l.h.s., so we obtain the growing-mode solution 5 m oc a in the deep matter era 
(H'/H ~ -3/2 and ~ 1). 

In GLPV theories, there is the anisotropic stress (77 7 1 ) between the gravitational potentials induced by the 
difference of c 2 from 1. Unlike the usual modified gravity models, the parameter 77 deviates from 1 even in the early 
cosmological epoch. We also note that the existence of the first term on the r.h.s. of Eq. (16.261) leads to the modified 
growth of matter perturbations relative to the case c 2 = 1 . 

In order to understand the evolution of the velocity potential, we take the AT derivative of Eq. (16.231) and then use 
other equations of motion. The resulting second-order equation for V m reads 


where 


ai = 


76 

4xi l 
2 


a 2 — (c t 1) 0 9 


Vm + a lVm + a 2V m — - ( 


4A(1 — f lm — x 2 ) — 76 xi (2 — f l m — 2x 2 ) 


Hn 


x> 


2x^ \aH 

3xi{4xf + 2(2 Q m — 3)x 2 -I- S2 m (fl m — 1)} — 76A(fl m -I- 4x 2 )(^m + x 2 — 1) 


2 xi L 

The general solution to Eq. (16.271) can be expressed in the following form 

V — V^ 

y m — r m -I- V m , 


(6.27) 


(6.28) 


(6.29) 


(6.30) 


where 177 is the special solution and Vm' 1 is the homogeneous solution derived by setting the r.h.s. of Eq. (16.271) to 
be 0. As long as the first term on the r.h.s. of Eq. (16.291) dominates over the other terms, the special solution is given 
by 


V w~-J_“kv = 

m — ‘2 A 


rx- 


r 2 _ IQ 
1 "t x “ L m 


(6.31) 
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In the second equality we used the fact that, for = 0, the sound speed squared c 2 is exactly given by Eq. (16.161) 
with u> = 0. 

The oscillations of perturbations are induced by the homogenous solution Vm \ In order to see the behavior of 
oscillations, we consider the scaling solution during the matter-dominated epoch. The scaling matter era can be 
realized by the dominance of the potential Uie _Al ^/ Mpl in Eq. (16.181) . which corresponds to jH] 

Xi=X 2 = ^ , Q m = l- (6-32) 

Substituting Eq. (16.321) into Eqs. (16.281) - (16.291) with A = Ai, we find that the homogenous solution satisfies 



(6.33) 


During the scaling matter era the scale factor evolves as a oc f 2 / 3 , so the evolution of the quantity 



(6.34) 


is known as K(Af) = Kje^/ 2 . where Ki is initial value of K at N 


0. Then Eq. (16.331) can be expressed as 


where 


U w ' 


c 2 j^-2 A/ - . 9(7A 2 -24) 
C °* Kie + 16A 2 


uL h) = o, 


U™ = a 3 / 4 V^ , 


C eff = 3 (A? — 3)(1 — c t) ■ 


(6.35) 


(6.36) 


If c 2 > 1, then we have c 2 ff < 0 and hence the perturbation Um' 1 is prone to the Laplacian instability on small scales. 
When c 2 < 1, the perturbation exhibits oscillations induced by the positive Laplacian term c 2 s K 2 (J\f ). Note that c 2 ff 
is related to c 2 as c 2 ff = c 2 — c 2 . As long as c 2 3> c 2 , c 2 ff is approximately equivalent to c 2 . 

The solution to Eq. (16.351) in the regime c 2 < 1 is given by 


= ciJ„(x) + c 2 Y u (x) , 


(6.37) 


where Ci,C 2 are integration constants, J„( x) and Y v {x) are the Bessel functions of first and second kinds respectively, 
with 


3a/7A? - 24 

v = —-- - - 1. 

2Ai 


X = 2 Ceff K (AT) = 


2c e ffk 

aH 


(6.38) 


Provided that x 1, the first term in the square bracket of Eq. (16.351) is much larger than the second term (which 
is of the order of 1). In this case the solution to Eq. (16.351) is given by Eq. (16.371) with the index v ~ 0. Using the 
asymptotic forms of the Bessel functions in the limit x 1, we obtain the following approximate solution 



(6.39) 


The perturbation Vm^ exhibits a damped oscillation with the frequency determined by c e gk. Since x oc e^/ 2 in the 
scaling matter era, the amplitude decreases as \Vm \ oc a -1 . As we will see later, the homogenous solution (16.311) 
stays nearly constant during the matter era. This means that, as we go back to the past, the oscillating solution 
(16.391) tends to dominate over the homogenous solution (as it happens for the dark energy models in f(R) theories 
0 )- The dominance of the oscillating solution can be avoided for the initial conditions satisfying {Vm 1 -C |V^|. 
This amounts to choosing the initial conditions of V m close to the value (16.311) and — 0. 
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C. Evolution of perturbations 


Numerically we integrate the perturbation equations of motion (16.20D - (16.23I) with (16.241) in order to find the precise 
evolution of V rn . S m as well as the gravitational potentials \F and $. 

First, we study the case in which the special solution to Eq. (16.271) dominates over the homogenous solution (16.391) . 
notwithstanding that this requires fine-tuning of the initial conditions for the perturbations. The initial conditions 
corresponding to such a case are given by V m ~ —K 2 x/ot 2 and V ^ — 0. From Eq. (16.231) the latter condition implies 
that SN ~ (H'/H)V m , in which case the r.h.s. of Eq. (16.201) is small such that £' ~ 0. As in the case of GR, we 
also employ the initial condition %' = 0 . Then, for given S m , the initial conditions of £, X) and 8 are known 
accordingly. 

In Fig. |T]we plot the evolution of V m , — *F, and $ versus the redshift z = ao/a — 1 (a o is the today’s value of a) 
for c 2 = 0.5, Ai = 10, A 2 = 0, and V 2 /V 1 = 10 -6 . During the matter-dominated epoch, the background solutions are 
in the scaling regime characterized by Eq. (16.321) . Numerically we find that the perturbations V m , x, and SN stay 
nearly constant by the end of the matter era for the initial conditions explained above. On using the approximate 
relations (16.311) as well as ~ 0 and x' — 0 with H’/H = —3/2, the evolution of perturbations during the scaling 
matter era can be estimated as 


V m ~ 


C - 


- - !!> 4,, 


c 2 K 2 - 3fl x 


SN ^ ^mSm 


5 c 2 — 2 c 2 

0C S ZC t Q r 

i Lrn.On 


2 c 2 K 2 - 3 fix ’ 

Jl „2 


2 c 2 (c 2 K 2 - 3fix) 

Then the gauge-invariant gravitational potentials satisfy 

3(l + c 2 -c 2 ) 


X - 


T ~ -- 


--IniSr: 


2(c 2 K 2 - 3fix) 
so that the anisotropic parameter is simply given by 

T] ~ 1 + 


$ ~ 


c 2 K 2 - 3 fix 

5Cs - 2 Ct(l + c 2 ~ c 2 t ) 
2 c 2 (c 2 I < 2 - 3 fix) 


flmSm • 


flmSm 


5(1 ~~ c t )( c s ~ c t ) 

3c t 2 (l + C 2 - C t 2 ) 


(6.40) 


(6.41) 


(6.42) 


where c 2 = + (A 2 - 3)(1 - c 2 t )/ 3. 

We confirmed that the analytic solutions (16.401) - (16.411) show good agreement with the numerical results shown in 
Fig. [Q In Fig. [^] we plot the evolution of 77 for three different values of c 2 smaller than 1. The anisotropic parameter 
is nearly constant during the scaling matter era. As estimated by Eq. (16.421) . the deviation of c 2 from 1 leads to the 
values of 7 / larger than 1. This is the observational signature of our model manifest in CMB temperature anisotropies 
[58| . From Fig. [I] we find that the two gravitational potentials —’F and $ start to decrease after the onset of the 
cosmic acceleration (z < 1). As we see in Fig. [ 2 j this leads to the variation of 77 , which signals the end of the scaling 
matter era. 

Substituting the solutions (16.401) into Eq. (16.261) . the perturbation S m during the scaling matter era obeys the 
following equation 



3 Geff j- ^ „ 

2 -^mOm — U , 


G eff 


A ~(c 2 + 1 — c 2 ) + 3flx „ 
c 2 K 2 - 3fl x 


As long as c 2 K 2 3> f2x, the effective gravitational coupling reduces to 


(6.43) 


Geff — 




(6.44) 


When c 2 < 1 we have G e ff > G. On using Eq. (16.441) . the growing-mode solution to Eq. (16.431) is given by 


Sm oc aP , 



' 1 + 24 


+ 24D m 



1 

4 ' 


(6.45) 


Compared to the case of GR (c 2 = 1), the growth rate of matter perturbations gets larger for c 2 < 1. However, since 
c 2 = c 2 + (Af — 3)(1 — c 2 )/3 and Ai > 10, this modification is suppressed to be small. For example, even for c 2 = 0, 
we have 24D m (l — c 2 )/c 2 = 72/A 2 < 0.7 in Eq. (16.451) . so the growth index 73 is not very different from that for c 2 = 1. 
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Figure 1 : The evolution of the velocity potential V m and the gravitational potentials — \P and $ versus 1 + z (= ao/a ) for 
c± = 0.5, Ai = 10, A 2 = 0, and V 2 /V 1 = ICE 6 . The present epoch (z = 0) is identified as f2 m = 0.3. We choose the initial 
conditions as Xl = x 2 = \/ 6 /( 2 Ai), x 3 = 0 , C = 1-4832 x 10 " 5 , x = -2.9124 x 10 “ 6 , 5 = 1.9351 x 10 " 3 , and V m = 1.8007 x 10" 7 
for the wave number K = k/(aH) = 25 at the redshift z = 527.9. These correspond to the initial conditions with the negligible 
oscillating mode characterized by V m = —K 2 x / a 2 and Vm = 0 . 


In Fig. [3] we show the evolution of fas for three different values of cf. When < 1 the effective gravitational 
coupling (16.44[) is bigger than G, so the growth rate of S m tends to be slightly larger than that for cj? = 1. However, 
the difference between the two cases cjr = 0.1 and cjf = 1 is very tiny except for the low-redshift regime in which the 
cosmic acceleration occurs. 

In the super-luminal regime c;( > 1, the Laplacian instability associated with negative can be avoided for 

c-t 1 < 77 ^ • (6-46) 

During the scaling matter era, this condition translates to Cj - 1 < 3/(Af — 3). The dashed curve shown in Fig. [3] 
(cj = 1.02) is the case in which the stability condition (16.461) is satisfied. In the low-redshift regime, fas for c% = 1.02 
is slightly smaller than that for c 1 = 1, so the gravitational force tends to be weaker. Numerically we confirmed that, 
for c 1 violating the condition (16.461) . the perturbations are prone to violent negative instabilities. 

Finally, we discuss the case in which the initial condition \Vm ‘ l | -C |I4n^| is violated. The numerical simulation 
shown in Fig. [I] corresponds to such an example, where the initial value of V m is not close to the special solution 
(16.311) . When cjf < 1, the velocity perturbation V m exhibits an oscillation induced by the sound speed c s during 
the scaling matter era. The amplitude of the homogenous solution decreases as \Vm \ oc a -1 , so V m approaches the 
special solution (16.311) in the end. Provided that the approach to occurs for z » 1, the evolution of fas in the 
low-redshift regime is similar to that shown in Fig. [3] 

In the numerical simulation of Fig. 01 we find that the gravitational potential — ’P shows a heavy oscillation around 
0 with a large amplitude. We recall that ’P is related to the perturbation SN = —V^ + (H 1 / H)V m , as T = 5N — 

[H'/H)x + x'■ From Eq. (16.391) the derivative Vm ^ can be estimated as 



+ C 2 cos ( x — 


(-;) 


(6.47) 


where we have assumed x = 2c e sK(J\f) 1. The amplitude of Vm' 1 is x/2 times as large as that of Vm' 1 ■ This is 
the reason why the amplitude of the oscillating mode of — ’P is much larger than that of V m . The definition of the 
anisotropic parameter 77 = — <P/\P loses its validity whenever — >P crosses 0, so in such cases we should directly resort 
to the effective gravitational potential (15.101) rather than rj. Furthermore, as — \P crosses 0, the effective gravitational 
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Figure 2: The evolution of the anisotropic parameter r/ versus 1 + 2 for cf = 0.1, 0.5, 0.9 with the model parameters Ai = 10, 
A 2 = 0, and V 2 /V 1 = 1CT 6 . The initial conditions are chosen to be similar to those given in Fig. [l] such that \Vm^\ < |V^|. 



1 + Z 


Figure 3: The evolution of fas versus 1 + 2 for cf = 0.1,1, 1.02 with the model parameters Ai = 10, A 2 = 0, and V2/V1 = 10 
The initial conditions are chosen to be similar to those given in Fig. [l] such that \Vm^\ -C ||. 


coupling G e g derived from Eq. (E3) changes its sign, leading, in this case, to a rather unclear interpretation for 
the gravitational interaction. In the numerical simulation of Fig. [4j the oscillation of —Hi does not damp away even 
around today (2 ~ 1). The amplitude of the oscillating mode of $ = ( + x is smaller than that of — T, such that the 
oscillation of the former disappears in the early stage of the matter era (see Fig. 0. 

The oscillation of the gravitational potential —HI leaves a distinctive imprint in CMB temperature anisotropies, so 
it should be possible to put tight constraints on the parameter space of initial conditions. We leave the analysis of 
observational constraints on the model (including the case of initial conditions | Vm <C |V^|) for a future work. 
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Figure 4: The evolution of the velocity potential V m and the gravitational potentials — U' and $ versus 1 + z for c? = 0.5, 
Ai = 10, A 2 = 0, and V 2 /V 1 = 1CF 6 with the initial value V m . = 5.0 x 1CF 7 and other initial conditions same as those in Fig. |T] 


VII. CONCLUSIONS 


In the framework of the EFT of modified gravity, we have studied the dynamics of cosmological perturbations and 
resulting observational consequences. Our starting point is the general action (13.11) in unitary gauge that involves a 
gravitational scalar degree of freedom and a matter field x mimicking a perfect fluid. This analysis encompasses the 
GLPV theories described by the Lagrangian (12.41) as a specific case. 

Expanding the action up to second order in the perturbations of ADM scalar quantities, we derived the scalar 
perturbation equations of motion (13.261) - (13.291) as well as the equation (13.361) for gravitational waves. In GLPV 
theories there is a non-trivial mixing between the propagation speeds of matter and gravitational scalar. The matter 
sound speed is not strongly affected by the deviation from Horndeski theories, but the modification to the scalar 
sound speed c s is large enough to be able to distinguish between GLPV and Horndeski theories. For non-relativistic 
matter with cj, = 0 we derived the exact relation — /3y, where ch is given by Eq. (14.71) and (3 h is a new term 

(14.91) arising in the theories beyond Horndeski. 

As an application of our general formalism, we studied the cosmology for a simple canonical scalar-field model 
described by the Lagrangian (EH) . The deviation from Horndeski theories can be quantified by the difference of the 
tensor propagation speed squared cjr = F((j>) from 1, where F(cj)) is a function of </> appearing in Eq. (16.21) . In the 
regime |q;h| = |c^ — 1| -C 1, the scalar sound speed squared is given by Eq. (16.161) . For dark energy models in 
which the ratio gets larger in the past (which is the case for most of models), |cj| can be much larger than 

1 even for |«h| "C 1. In particular, for cjf > 1, the scalar perturbation is plagued by the Laplacian instability unless 
the condition (16.171) is satisfied. 

The behavior of large values of in the past can be avoided for the dark energy model given by the potential 
(16.181) . in which case there exists the scaling solution characterized by the constant ratio H m /f lx = 2(A^ — 3)/3. In 
this model we studied the evolution of cosmological perturbations to find observational signatures in the framework 
of GLPV theories. The solution to Eq. (16.271) for the velocity potential V m = F[v m can be written in terms of the 
sum of the special solution Vm ' 1 and the homogenous solution Vm 1 . During the scaling matter era Vm^ stays nearly 
constant, whereas Vm ' 1 decreases in proportion to a -1 with oscillations. 

Provided that |I4n^| |I4n^|, the evolution of perturbations during the scaling matter era can be known analyti¬ 

cally. Of course this requires the fine-tuning of initial conditions. In this case, the anisotropic parameter 77 = — <1>/’I' 
between the two gravitational potentials is expressed in terms of and cl? alone. As we see in Fig. [2j the deviation 
of 77 from 1 tends to be larger as cj? is away from 1. We have also estimated the growth rate of matter perturbations 
S m and found that the evolution of 5 m is not sensitive to the tensor propagation speed ranging in the regime (16.461) . 

If Vm' 1 dominates over Vm^ at the initial stage of the matter era, the perturbations exhibit rapid oscillations with 
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the frequency related to c 2 until V m approaches the special solution Vm* ■ The amplitude of oscillations is particularly 
large for the gravitational potential T, see Fig.Q] This property conies from the fact that 'F is related to the derivative 

term Vm* , whose amplitude is x/2 times as large as that of Vm*. 

We have thus shown that the simple extension of Horndeski theories to GLPV theories gives rise to interesting 
observational signatures. For the initial conditions satisfying |V^| <C |V^|, the anisotropic parameter tj deviates 
from 1 even in the early stage of the matter era (which is usually not the case in Horndeski theories). For other initial 
conditions, the gravitational potential 'F rapidly oscillates with a large amplitude especially when c 2 is away from 1 . 
In both cases, we should be able to put tight bounds on the deviation parameter an = c 2 — 1 from the observations 
of CMB and weak lensing. It will also of interest to study local gravity constraints on our model along the lines of 
Refs, dm 
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Appendix A: Quasi-static approximation on sub-horizon scales in Horndeski theories 


In this Appendix we derive the expression of G e ff and <F e g in Horndeski theories under the quasi-static approximation 
for the modes deep inside the sound horizon. Since an = 0 in this case, we can solve Eqs. (15.1211 and (15.1411 for if) and 
C, as 


where 


£,s(46! - c t 2 W) hW-hLg 

hctL^s - 4 b\ ’ b 2 (%L' S ~ 46 ? 


6 i = L iS + HL tS , b 2 = W + HW + p. 


(Al) 


(A2) 


Substituting the relations (I A 111 into Eq. (15.1311 with the approximation S m — <5, we obtain Eq. (15.711 with the effective 
gravitational coupling 


G eff 


_fr 2 c 2 £,5 ~ 4 b\ _ 

4tt£ s (W 2 c 2 + 46 2 £,5 - 8&i W) ‘ 


(A3) 


From Eq. ED the gravitational potential $ = £ + Hip satisfies the relation $ 


— 77 T, with the anisotropic parameter 


_ 6 2 £,g - 61W + HL t s(Wc$ ~ 46 i) 

77 b 2 c^L $ - 4b\ 

On using Eq. (ism . ED. and (IA4I) . the effective gravitational potential (15.101) obeys 


(A4) 


a 2 6 2 (1 + c 2 )£,5 — 4 b\ - biW + F7£,s(Wc 2 — 4&i) 
fc 2 2£, s (W 2 c 2 + 46 2 £, 5 - 861 W) P m 


(A5) 


The results (IA3I) - (IA5I) match with those derived in Ref. 0 in the Newtonian gauge by taking the massless limit 
m<j, —> 0 . 
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